Plume motion and large-scale circulation in a cylindrical Rayleigh-Benard cell 
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We used the time correlation of shadowgraph images to determine the angle O of the horizontal 
component of the plume velocity above (below) the center of the bottom (top) plate of a cylindrical 
Rayleigh-Benard cell of aspect ratio F = D/L — 1 (D is the diameter and L ~ 87 mm the height) 
in the Rayleigh-number range 7x10' < B < 3 x 10 9 for a Prandtl number a = 6. We expect that 
O gives the direction of the large-scale circulation. It oscillates time-periodically. Near the top and 
bottom plates Q(t) has the same frequency but is anti-correlated. 

PACS numbers: 47.27.-i, 44.25.+f,47.27.Tc 



Turbulent Rayleigh-Benard convection (RBC) is an 
important process that occurs in the oceans, the atmo- 
sphere, the outer layer of the sun,the Earth's mantle, 
and in many industrial processes, [lj Despite the seem- 
ing simplicity of the idealized laboratory experiment, 
i.e. a fluid between horizontal parallel plates heated 
from below and cooled from above, our understanding 
of the physical mechanism of turbulent RBC remains 
incomplete. 0, 0, 01 Much of the heat transport is me- 
diated through the emission of hot (cold) plumes of fluid 
from a thin boundary layer adjacent to the bottom (top) 
plate, and these plumes are carried by a large-scale circu- 
lation HI1HB11HII3 known as the " wind of turbu " 
lence". For systems with aspect ratio T = D / L = 0(1) 
{D is the diameter and L the height of a cylindrical cell) 
this wind, when time averaged, takes the form of a sin- 
gle convection roll filling the entire cell. The plume in- 
teraction with this circulation is a central component of 
the turbulent RBC problem, and yet only little is known 
quantitatively about plume motion and the wind. To a 
large extent we expect that at least the horizontal mo- 
tion of the plumes is slaved to the wind velocity, since 
the only independent force acting on the plumes is the 
buoyancy force in the vertical direction. Here we present 
a quantitative study of plume motion and thus, by in- 
ference, of the wind direction above (below) the center 
of the bottom (top) plate. The measured horizontal di- 
rection oscillates time periodically, and the oscillations 
persist with a unique frequency over hundreds of cycles. 
The oscillations at the top and bottom have the same 
frequency but a phase which is displaced by half a cycle. 
We conclude that the wind is a significantly more compli- 
cated dynamical system than a simple Rayleigh-Benard 
convection cell. 

Numerous measurements of the speed of the fluid 
and/or of the temperature at distinct points (i.e. of 
scalar quantities) were made before and revealed a time- 
periodic component. 0, HI El El El I n some instances 
these results were interpreted as indicative of a time- 
periodic plume emission from the plates. Where direct 
comparison is possible, we find that the measured fre- 



quencies agree quantitatively with our determination of 
the oscillation frequency of the horizontal wind direction. 
It has been shown 9] that the corresponding oscillation 
period is commensurate with the period of circulation 
of the wind. It had been suggested that periodic plume 
emission from the top and bottom plates, synchronous 
with the wind circulation period, provided a mechanism 
for driving the wind. Our results offer a different expla- 
nation of the periodicity of the system, and thus there is 
no longer a direct reason to assume that periodic plume 
emission is a central ingredient of the driving mechanism. 

Rayleigh-Benard convection is characterized by the 
Rayleigh number R = agATL 3 /Ki/ and the Prandtl 
number a = v j K. Here a is the isobaric thermal ex- 
pansion coefficient, k the thermal diffusivity, v the kine- 
matic viscosity, g the acceleration of gravity, and AT the 
applied temperature difference. We used methanol with 
a = 6.0. For 0.5 < AT < 20°C at a mean temperature of 
40°C they covered the interval 7 x 10 7 < R < 3 x 10 9 , in a 
range sometimes called the " hard turbulence" regime [|| , 
but below the Kraichnan regime [l^ where the boundary 
layers themselves become turbulent and break up. 

The apparatus was described before. 

III El The 

cell had a sapphire top plate of thickness 0.32 cm. The 
bottom plate, of aluminum 0.65 cm thick, had a top sur- 
face which was polished to a mirror finish and was heated 
from below by a metal-film heater uniformly distributed 
over the cross section of the cell. The wall was made 
of high-density polyethylene. The internal cell diameter 
and height were D ~ 87.4 and L ~ 86.6 mm respectively. 

Above the cell was a shadowgraph optical system El • 
Shadowgraph images represent a vertical average of the 
horizontal Laplacian of the refractive index (and thus 
temperature) field in the cell. In the central part of 
the cell, between the thermal boundary layers, the tem- 
perature undergoes vigorous fluctuations but on average 
is nearly equal to the mean temperature of the system. 
This part of the cell does not contribute significantly to 
the shadowgraph image. The laminar boundary layers 
themselves also do not contribute because they have no 
significant horizontal temperature variation. Thus the 



2 




17.2s 



FIG. 1: Shadowgraph images of plumes for R = 5.6 x 10 . 
The elongated black stripes are the hot plumes just above the 
bottom plate. The white squares in the images, corresponding 
to 31 x 31mm 2 , indicate the part used to calculate the tem- 
poral correlation between two successive images separated in 
time by 1.22s. The auto-correlation functions are shown on 
the right. A black dot at their center marks their origin. 



major contribution to the images comes from the plumes 
which emerge from the bottom or top boundary layers. 
These consist of localized relatively warm or cold fluid 
respectively, and thus appear as dark or bright areas. In 
order to reduce contributions from optical imperfections 
of the system, the images are divided by a reference im- 
age taken with AT = and then rescaled. Two typical 
divided images are shown in the left column of Fig. ^ 
With the particular setting of the optics, the warm dark 
plumes at the bottom are more apparent than the cold 
bright ones at the top, but both are present. In an at- 
tempt to break the rotational invariance of the system 
and to give the large-scale flow a preferred direction, our 
measurements initially were made with the entire system 
tilted by an angle of 2°. However, a subsequent sequence 
of measurements with the system horizontal within 0.1° 
gave the same results. Using camera exposure times of 30 
or 60 ms, we took sequences consisting of 8192 or 16386 
images at time intervals of 200 to 500 ms, correspond- 
ing to an acquisition time of 30 to 150 min or 50 to 350 
periods of the oscillations, for a given R. 

Plumes were visualized from the side m any times, ei- 

" J! IU or by 
Those ob- 



ther by the shadowgraph method ^ 




si 



using thcrmochromic liquid crystals 
servations suggest that plumes have the shape of a mush- 
room with a cap, and with a stem connecting the cap to 
the boundary layer. We find that visualization from the 
top reveals new features. 01 -"- n Fig. ^ one sees that the 
plumes are extended horizontally, consistent with rising 
or sinking sheets of fluid. We only very rarely observed 
a typical mushroom cap at the lateral termination of a 
plume. The plumes in our experiment have a horizontal 
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FIG. 2: Temporal evolution of the angle O of the horizontal 
velocity component of the plumes close to the bottom plate. 
As for Fig. |TJ R = 5.6 x 10 8 . 



extent that varies from a few mm to a few cm, and a 
thickness of around 1 to 2 mm for R = 2 x 10 8 . This 
plume thickness is about the same size as the thermal 
boundary-layer thickness, as expected. jU Compar- 
ison of successive images, separated in time by a second 
or so, shows that the plumes are moving, somewhat errat- 
ically, but on average in the direction of their long axis. 
The alignment of convection rolls by a prevailing wind or 
by shear flow is a well known phenomenon. |2 ll |25j . 

Visual observation revealed that the direction of the 
plume movement oscillated in time. We calculated 
time auto-correlation functions between successive im- 
ages. First, a threshold was applied and a binary image 
was obtained, rendering the dark plumes in black. Then, 
on a central part of the image of size 31 x 31mm 2 , illus- 
trated by the squares in the images of Fig. ^ the cor- 
relation between successive images was computed. Two 
examples of these correlation functions are shown in the 
right column of Fig. ^ The vector from the origin to 
the highest peak gave the average displacement of the 
plumes which occurred during the time interval between 
the 2 images, and thus the plume velocity. Here we focus 
on the angle of this vector relative to an axis which 
was chosen so that the time average of vanished. 

Results for are shown in Fig. [21 As already noted 
by visual observation, oscillated in time with a well 
defined frequency to. The oscillations covered an angular 
range of more than 90°. Results for u>, scaled by the verti- 
cal thermal diffusion time L 2 /n = 7.88 x 10 4 s, are shown 
in Fig. as a function of R for both the slightly tilted 
(open circles) and the untilted (solid circles) apparatus. 
The two sets agree quite well with each other, and can 
be described by ui = uioR -13 , with luq = 1.18(1.26) and 
(3 = 0.465 ±0.011(0.460 ±0.012) for the tilted (untilted) 
run. 

The existence of time-periodic behavior in our range 
of R was observed previously, by measurements either of 
the temperature or of a component of the velocity (i.e. 
of scalar quantities) at isolated poin ts in the convection 
cell, in cylindrical |6J, |8J, |9j, |lfj, |llj as well as rectangu- 
lar |l9| systems. These measurements could determine 
the frequency, but could not identify the nature of the 
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FIG. 3: The oscillation frequency to, in units of the vertical 
thermal diffusion time L 2 /k, of the angle Q of the plume 
motion just above the bottom plate for the slightly tilted 
(open circles) and the untilted (solid circles) apparatus. The 
straight line correspond to a fit of a powerlaw to the data of 
Ref. 0- 



oscillating mode. The measured frequencies depended 
on R (as shown in Fig. [3J and a. For the same a, our 
results agree quantitatively with those obtained by the 
scalar measurements, 9] as illustrated by the solid line in 
Fig. |2l which represents alitw = 1.26i? ' 46 to data with 
a = 5.4. J9| Thus the scalar and our vector measurements 
capture the same physical phenomena. In some cases 
this frequency was attributed to periodic plume emission 
[ill Eil Ejf . Our results indicate that its origin is the 
horizontal oscillation of the large-scale flow. 

It had been determined by others that, consistent 
with a model given in Ref. |27| . the frequency is pro- 
portional to the Reynolds number R e of the large-scale 
circulation. H H ES 0H The prediction for R e by 
Grossmann and Lohse |3l| . for our a and over our range 
of R, gives (3 as an effective exponent with a value of 
0.446, slightly lower than but reasonably consistent with 
our measurements and with those of Ref. || . 

We obtained the auto-correlation function 



C„(t) = (5G(t)6e(t + T))/(5Q(t) 2 ) 



(1) 



from image sequences covering typically 8000s. Here 
SQ(t) = Q(t) — (6(i)) (angular brackets indicate a time 
average). An example in Fig. 01 shows some decay for 
the first 500s, but then continues to oscillate indefinitely 
with an amplitude which varies irregularly. This indi- 
cates that the flow maintains its characteristic frequency 
indefinitely with high accuracy, but does not rule out a 
slow phase shift with time. 

The shadowgraph apparatus had a depth of field large 
enough to visualized the hot plumes near the bottom 
(which appear dark), and the cold plumes near the top 
(which appear bright). By choosing different threshold 
values in the image processing, we were able to sepa- 
rate the hot and cold plumes, and determine the velocity 
of both of them from the same images. Both velocities 
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FIG. 4: Auto-correlation of the angle 9 for R = 8.37 x 10 8 . 
The lower panel shows a small time interval of the upper one 
on an expanded scale. 
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FIG. 5: Cross correlation between the direction of the hot 
and cold plumes for R = 9.77 x 10 s . 



have similar oscillation amplitudes and the same lj. An 
example of the cross-correlation function of the hot- and 
cold-plume angles Oh and C 



C c (r) 



(8Q h (t)5e c (t + T)) 

y/(s&i(t))(se 2 c(t)) 



(2) 



is given in Fig. [S] and shows that the directions of motion 
of the hot and cold plumes are anti-correlated, i.e. out 
of phase by ir. This was the case for all R. 

Let us summarize our understanding of the large- 
scale flow. One of the interesting features of Rayleigh- 
Benard turbulence began to emerge when it was found 
that the flow, although vigorously fluctuating, highly ir- 
regular, and seemingly random, contains a well defined 
frequency. This frequency is Rayleigh-number depen- 
dent, and its period is proportional to a turnover time 
of a large-scale circulation. The Reynolds numbers of 
the large-scale flow implied by the frequency measure- 
ments agree well with the prediction of Grossmann and 
Lohse. [3l| The existence of the oscillations of the flow 
direction shows that the flow is much more structured 
than initially assumed. Due to contributions of a num- 
ber of investigators, HSU H El El El SH3 a lot of 
progress was made in our understanding of the large-scale 
circulation, its periodicity, the boundary layers, plume 
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FIG. 6: Schematic summary of a snapshot of the plume move- 
ment. The small angle at the bottom labeled 2° indicates the 
tilt of the apparatus relative to gravity. The curved double- 
headed nearly horizontal lines near the top and bottom indi- 
cate the amplitude of oscillation of the angle 0. The arrows 
through the top and bottom plane are an example of the di- 
rection of flow at an angle <d at the top and — O at the bottom 
at a given time. Although not directly evident from our mea- 
surements, we assume that the instantaneous direction of flow 
near the side wall is given by connecting the head of the top 
(bottom) arrow with the tail of the bottom (top) arrow. 



dynamics, and the plume interaction with the large-scale 
wind. A schematic drawing of the velocity field and the 
plume distribution had evolved Q and was named the 
"flywheel" by some. It showed the fluid circulating as a 
loop in a well defined plane with a vertical axis, station- 
ary in time, and the plumes rising and falling on opposite 
sides of the loop. Alhough this picture was appealing, it 
was called into question already by measurements of the 
characteristic frequency and of the amplitude of the os- 
cillations at various places in a cell.|26j 

Our work focused on the plumes viewed from the top. 
Our first result is that near the bottom and the top plate, 
plumes when viewed from above look like sheets. They 
move across the field of view, and are elongated in the di- 
rection of their movement. Our second result is that the 
frequency measured in the temperature or flow field cor- 
responds to a periodic oscillation in the horizontal plane 
of the direction of the plume movement, in contradic- 
tion to the simple flywheel model. Our third result is 
that the direction of the hot and cold plumes near the 
top and bottom plates are out of phase by tt. A sim- 
ple sketch summarizing the new features of the flow is 
given in Fig and shows that the schematic picture of a 
simple loop (flywheel), partly driving and in turn partly 
driven by plumes, is an over-simplification. It turns out 
that the large-scale flow has a more complicated three- 
dimensional spatial structure, and that it may be re- 
garded as an interesting dynamical system with a char- 
acteristic frequency of a three-dimensional mode. 

This work was supported by the US Department of 
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